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ABSTRACT 

In galaxy clusters the entropy distribution of the IntraCluster Plasma modulates the latter's equilibrium within the Dark Matter gravi- 
tational wells, as rendered by our Supermodel. We argue the entropy production at the boundary shocks to be reduced or terminated 
as the accretion rates of DM and intergalactic gas peter out; this behavior is enforced by the slowdown in the outskirt development at 
late times, when the Dark Energy dominates the cosmology while the outer wings of the initial perturbation drive the growth. In such 
conditions, we predict the ICP temperature profiles to steepen into the cluster outskirts. The detailed expectations from our simple 
formalism agree with the X-ray data concerning five clusters whose temperature profiles have been recently measured out to the 
virial radius. We predict steep temperature declines to prevail in clusters at low z, tempered only by rich environs including adjacent 
filamentary structures. 

Key words, galaxies: clusters: general — galaxies: clusters: individual (Abell 1795, PKS 0745-191, Abell 2204, Abell 1413, Abell 
1689) — X-rays: galaxies: clusters — methods: analytical. 



1. Introduction 

Galaxy clusters constitute the largest bound structures in the 
Universe, with their masses up to M ~ IO'^Mq and outskirts 
extending out to sizes R ~ a few Mpcs. These set the interface 
• between the intergalactic environment keyed to the cosmology 
I at large, and the confined intracluster plasma (ICP). The latter 
. pervades the clusters at temperatures kgT oc GM/R ~ 5 keV and 
' number densities n ~ 10"^ cm"-', and so emits copious X-ray 
powers mainly by thermal Bremsstrahlung (see Sarazin 1988). 
The ICP coexists with the gravitationally dominant dark mat- 
' ter (DM) component in the baryonic fraction m/M close to the 
cosmic value 0.16, and the two build up together from accretion 
across the cluster boundary. 

The build up comprises an early collapse of the cluster 
body, tailing off into a secular developmen t of the outskirts 
, by smooth accretion and minor mergers (see Zhao et al] 20031: 
'. iDiemand et ani2007t IVass et al.ll2008t iNavarro et al.ll2010l) . In 
' radius, the body ranges out to r ~ r_2 where the slope of the DM 
density run n(r) equals -2; the adjoining outskirts extend out to 
the current virial radius R with steepening density. 

In time, the transition is marked by the redshift z,; there- 
after r_2 stays put while R grows larger in a quasi-static, 
self-gravitating DM equilibrium (described through the Jeans 
equation, see Lapi & Cavaliere 2009a), to imply for the stan- 
dard concentration parameter c = R/r-2 observed values 
c ^ 3.5 H(zi)/H(zobs) in terms of the Hubble parameter H{z). 
In the following w e adopt the standard flat cosmology (see 
iDunklev et al.ll2009l) . So values of c ranging from 3 to 10 cor- 
respond for Zobs ~ - 0.2 to young or to old dynamical cluster 
ages Zi ~ 0.2 - 3. The concentration can be directly if labo- 
riously probed with gravitational lensing (see iBroadhurst et al.l 
[2008; Lapi & Cavaliere 2009b|). 

Secular accretion of DM goes along with inflow of inter- 
galactic gas. The ensuing ICP equilibrium is amenable to the 



powerful yet simple description provided by the SupermodeQ 
(SM; see Cavahere et al. 2009, hereafter CLFF09). 

Clearly, inflows into the cluster outskirts are exposed to the 
cosmological grip. This is the focus of the present paper 

2. Entropy run vs. cluster buildup 

The SM expresses in full the hydrostatic equilibrium (HE) of the 
ICP in terms of DM gravity and of the 'entropy' k = ksT/n^^^. 
In its basic form, the latter's physical run may be represented as 



k(r) = k, + {kK-k,)(r/Rr 



(1) 



consisten t out to r a; R/2 with recent analyses of wide cluster 
samples (ICavagnoIo et al. I l2009t iPratt et al.ll20ra) . This embod- 
ies two specific ICP parameters: the central level kc and the outer 
powerlaw slope a. 

The former is set at a basal level kc ~ 10 keV cm^ 
by in termittent entropy inj e ctions by central AGN feedback 



ICavahere et al.' '2002^ 

120001; pcNamai-a & Nulsenl 
condition corresponds to cool 



Valageas & SilS [19991 IWu et al.l 
2007h ; the ensuing quasi-stable 
morphologies (CC; 



core morptiologies (CC; see 
Molendi & Pizzolato 2001.) , featuring a limited central temper- 
ature dip and generally large concentrations c « 6 - 10, as dis- 
cussed by CLFF09. 

On the other hand, k^ may be enhanc ed up to sev- 
eral 10^ keV cm^ by deep merg ers (e.g., Mc Carthy et al.l 
120071: iMarkevitch & 'Vikhlininll200'7l) . frequent during ±e clus- 
ter youth; these events give rise to non cool core (NCC) clus- 
ters, featuring generally low concentrations c 3 - 5 and 
a central temperature plateau, scarred in some instances by 
imprints from recently stalled blast waves (see discu ssions by 
iFusco-Femiano et al.|[2009: .Rossetti & Molendill2010l) . 



' IDL and FORTRAN algorithms which implement the Supermodel 
and run in a fraction of a second on a standard laptop, can be found at 
http://people.sissa.it/~Iapi/Supermodel/. 
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Fig. 1. Examples of entropy profiles (normalized at r = r/j) and 
slopes (inset). Dotted line: after Eq. (1) with A:^ = and a - 1.3. 
Solid line: Eq. (4) with /t,. = 0, a = 1.3, = 0.25 7?, and a' =0.5. 
Dashed line: Eq. (4) with = 0, a = 1.3, = 0.25/?, and 
a' = 1.5. 



2. 7 . The outer regions 

The second term in Eq. (1) describes the powerlaw outward rise 
expected from the scale-free stratification of the entropy con- 
tinuously produced by the boundary accretion shock, while the 
cluster grows larger by slow accretion. 

The slope at r » /? with standard values around 1 has been 
derived by CLFF09 from the shock jumps and the adjoining HE 
maintained by thermal pressure, to read 

ur = 2.37 - 0.47 bu . (2) 

Here bji = fimpV^^/kBTu mar ks the ratio of the potential to the 
thermal energy of the ICP (see lLapi et alj2005l:IWtll2005h . This 
reads bR =s 3 v|/2 AO when a strong shock efficiently thermal- 
izes into 3 degrees of freedom the infall energy; the latter is ex- 

^ dr G 6M/r , 

experienced by successive shells of DM and gas that expand; 
owing to the excess mass 6M they turn around at the radius 
7?ta ~ 2 7? to start their infall toward the shock at R. 

At any given cosmic time t, Eq. (2) holds at the current 
virial radius R{t). On the other hand, the entropy deposited there 
is conserved during subsequent compressions of the accreted 
plasma into the DM gravitational well, while no other major 
sources or sinks of entropy occur down to the central 10^ kpc. 
Thus while the cluster outskirts develop out to the current ra- 
dius R{t), the specific entropy stratifies with a running slope 
a(r) = aR(t) that retains th e sequence of original values set at 
the times of deposition (see lTozzi & Normanll2001h . 

Values fl X! 1 obtain on adopting the standard ratio R/Rt^ 
^ 0.5, and the simple potential drop AO/v^ ^ 1 - (R/Rta) ~ 0.5 
associated to a flat initial mass perturbation 6M/M cc M"*^ with 
e ~ 1 that describes the collapse of the cluster body as a whole. 
This implies bR ^ 3 and a ^ I from Eq. (2). In fact, AO/v^ x 
0.57 obtains from the full DM a-profile (see CLFF09), implying 
bR ~ 2.7 and the standard value a » 1.1. 

However, considerable variations of a{r) are to be expected 
in the cluster outskirts, as discussed below. 



2.2. Development of the outskirts 

The cluster outskirts for r > r_2 originate from the wings of a 
realistically bell-shaped perturbation, that may be described by 
SM/M oc M""^ for e exceeding 1 (e.g.. .Lu et al..,2006.) . Then the 
outer potential drop 

AO) 1 - (R/RJ''-^ 



is shallower relative to the body value, so leading to higher val- 
ues of bR and lower values of a. Thus as r increases outwards we 
expect k(r) to deviate downward from a simple powerlaw. 

The argument may be phrased in terms of the accretion rate 
M. A shell 6M enclosing the mass M will collapse when 6M/M 
attains the critical thres hold 1.686 £)~'(f ) in terms of the linear 
growth factor D(t) (e.g.. I Weinbergll2008h . So the shape parame- 
ter e also governs the mass buildup after M oc D'^'^ oc f''^^; here 
we have represented the growth factor as D(t) oc f' with d rang- 
ing from 2/3 for z ^ 1 to approach 1 /2 as z ^ 0. So the outskirts 
develop from the inside-out, at accretion rates M/M d/et 
that lower for e exceeding 1, and for d decreasing toward 1/2 
at late cosmic timefl We add that at small ac cretion rates the 
shock position outgrows R (seeVoit et^alj2003|), while the shock 
strength may weaken; both these effects will decrease a relative 
to Eqs. (2) and (3), and will be discussed in detail elsewhere. 

So we see that flatter slopes a of the entropy are to prevail 
for decreasing accretion rates M of DM and gas; these have a 
twofold origin. First, the cosmological structure growth slows 
down at later cosmic times (low Zobs). as expressed hy d < 2/3. 
Second, perturbation wings marked by e > 1 imply shallow 
gravitational wells and little available mass to accrete in average 
environs; the effect may be locally offset (and represented with a 
smaller effective e) in specifically rich environments, including 
adjacent filamentary large-scale structures. 

The decline of a from the body value and the entropy bend- 
ing set in at a radius r_2 where matter began stratifying onto 
the outskirts just after z,. Such a radius is evaluated in terms of 
the observed concentration in the form rb/R ~ r-2/R ~ 1/c, 
to take on values around 0.2 - 0.3 for typical concentrations 
c X 6 - 8 of CC clusters. H ints to this trend loo m out in the 
data bv Pratt et all (1201 Ol) and lHoshino et al.1 (1201 Oh . 

To sum up, under the lower accretion rates prevailing at later 
times in average environs, we expect the entropy run to flatten 
out or even decline into the cluster outskirts; then the tempera- 
ture will decline as kgTir) = kn^^^ oc r^^ or steeper, after Eq. (7) 
of CLFF09. Do such behaviors show up in real clusters? 

3. A case study on current data 

Toward an answer, we use the SM to provide profiles of density 
and temperature from expressing the expected entropy run in a 
simple form: the initial slope a still applies in the cluster body 
for r < rt,; but for r > it goes over to a decline toward the 
current boundary value aR < a following a - a' {r j rt, - 1) with a 
constant gradient a' . So the entropy profile reads 

i kc + (kt - kc) (r/rbY r < rb 
k(r) = i . (4) 

[kRir/Ry^"' e"'^'^-'''''''' r>rb 

- Note that d — » would occur if the Dark Energy de nsity increased 
with time to cause an ultimate Cosmic Doomsday (see ICaldwell et all 
2003); this would imply truly vanishing accretion rates, and result into 
a cutoff of the temperature profiles. 
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Abell 2204 @ z = 0.15, R. 




Fig. 3. Same as Fig. 2 for the low-z clusters PKS 0745-191 (top) 
and A1795 (bottom). Data are from lSnowden et al. (2008) with 
XMM-Newton (blue circles) and from lGeorge et al.l (l2009() and 
Bautz et al. (2009) with Suzaku (red squares). 

Table 1. Fitting parameters from ICP temperature profiles 
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Fig. 2. Profiles of projected X-ray temperature (brightness in • DM concentration estimated from X rays. 

the insets) for the CC clusters A1689 (top), A1413 (mid- " Reduced values for fits with the entropy run in Eq. (4); in paren- 

dle) and A2204 (bottom). Data are from Snowden et al. ('2008') thesis the values with a' = 0, corresponding to Eq. (1)_^ 

with XMM-Newton (blue cir cles) and fro m Kawaharada et al. ' CC classification is controversial fsee .Riemer-Serensen et alJ200^ . 
(120 Id) . iHoshino et all (l2010l) . and iReiprich et all (2009) with 
Suzaku (red squares). Our best-fits with the SM after Eq. (4) 

are illustrated by the solid lines, while dashed lines refer to the 3 
fits based on Eq. (1). 
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The outer branch describes a linear decline of the slope with the 
gradient a' = (a - ai{)/(R/rb - 1); normalizations have been set 
as to obtain for k(r) a continuous function and derivative. We 
illustrate in Fig. 1 examples of entropy profiles after Eqs. (1) 
and (4), with parameters indicated by the analyses below. 

Adding to the several CCs and NCCs previously analyzed 
out to R/2 in terms of the standard entropy run (see CLFF09 and 
[pusco-Femiano et al. 2009), here we focus on five CCs with data 
now available out to r ^ R, and analyze them in terms of the 
entropy run of Eq. (4). We report in Table 1 the resulting param- 
eters with 68%-level uncertainties. For these old CC clusters we 
take kc ~ 10 keV cm^ as anticipated in § 2. 

4. Discussion and conclusions 

These data clearly bear out our expectations of steep runs for 
T(r), stemming from outer entropy production being reduced as 
the accretion rates M of DM and intergalactic gas peter out (see 
§ 2). We trace back such a behavior to two concurring sources: 
(i) the cosmological slowdown in the growth of outskirts devel- 
oping at late cosmic times; (ii) shallow perturbation wings scant- 
ily feeding the outskirts growth in average or poor environs. 

Under such conditions, we show that on average the entropy 
profiles are to flatten out or even decline into the cluster out- 
skirts as represented by Eq. (4). Then our Supermodel yields 
steep outer profiles of projected temperatures (and flatter densi- 
ties), in close agreement with current data out to R. 

We also expect the cosmological decrease of the accretion 
rates to be locally offset in rich environs biased high, or in clus- 
ter sectors adjacent to filamentary large-scale structures; more 
standard runs of temperature are pre dicted there. These loom o ut 
in a sector of A1689 as observed bv lKawaharada et aP (l2010h . 

Compared with numerical simulations like the ones by 

iNag ai et al. (2007) and Roncarelli et al. (2006), our temperature 
declines and entropy dearths are comparable as for clusters like 
A2204 and A1413; but they are quite sharper for low z clusters 
like A1795 and PKS0745-191, well beyo nd possible contami - 
nations due to Suzaku PSF smearing (see Rei prich et al.ll200"9l) . 
Low-noise, high-resolution simulations addressing the issue of 
entropy production related to richness of surrounding environ- 
ment or adjoining filamentary structures will be most fruitful. 

We have checked that with the SM the gas mass m(< r) 
grows monotonically outwards, and actually faster than the 
DM's, so as to produce an i ncreasing bar yonic fraction by a fac- 
tor 10 from R/20 to R (cf. IZhang et alfe oiO). So sharply de- 
clining T{r) and increasing masses are in fact consistent with 
thermal HE. Note that as M decreases to the point that the infall 
velocities decrease to transonic values, the shocks weaken, ther- 
malization becomes inefficient (see CLFF09), the entropy pro- 
duction is terminated, and thermal pressure alone cannot support 
HE any longer. 

Equilibrium in the outskirts can be helped by bulk or turbu- 
lent, merger-induced motions; these contribute up to 10 - 20% 
of the total support in relaxed clusters, as gaug ed in terms of 
pressure and X-ray masses both observationally (iMahdavi et al.l 
12008 ': 'Zhang et al.' '2010', adding to the refs. in § 3) and 
nume rically (Nagai et al. 2007; Piffaretti & Valdarnini 2008; 
iBode et all l2009t"lLau et al.1 120 lOt iMeneg hetti et alj l2010l) . A 
similar accuracy is intended for our predictions based on thermal 
HE, while a finer approximation with non-thermal contributions 
included in the SM will be provided elsewhere. 

On the other hand, even extended HE constitutes just a useful 
approximation, due to fail even in relaxed clusters: locally, when 
minor lumps of cold gas fall into the cluster from an adjacent 



filament, along with smooth accretion; globally, beyond^ where 
the DM equilibrium is in jeopardy. 

In sum, from the SM with thermal HE and optimal energy 
conversion in strong shocks, we find temperature profiles declin- 
ing sharply outwards. These stem from progressive exhaustion 
of mass inflow, especially for clusters in average or poor envi- 
ronments at low Zobs\ such a predicted trend is consistent with 
the current evidence for the handful of clusters in Figs. 2 and 3. 

Our punchline is that cosmology - besides aff'ecting the 
clus ter statistics like mass or temperature distributions (see 
Vik hlinin et al. I I2009I) - concurs with the perturbation shape to 
set the outer structure of individual clusters and their develop- 
ment. 

All these rich phenomena expected at the interface be- 
tween the ICP and the intergalactic medium call for ex- 
tensive probing even at z S 0.2 with the next genera- 
tion of X-ray telescopes planned to detec t at high resolutions 
low-surface brightness plasma (see iGiacconi et al. .2009 ), like 
WFXT (see http://wfxt.pha.jhu.edu/) and eventually IXO (see 
http://ixo.gsfc.nasa.gov/). 
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